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Abstract 



We consider the motion of a particle in a periodic two dimensional flow per- 
turbed by small (molecular) diffusion. The flow is generated by a divergence free 
zero mean vector field. The long time behavior corresponds to the behavior of the 
homogenized process - that is diffusion process with the constant diffusion matrix 
J> ■ (effective diffusivity) . We obtain the asymptotics of the effective diffusivity when 

, the molecular diffusion tends to zero. 

m 
q 

cn 1 Introduction 

o 

o\ . 

O . Consider the following stochastic differential equation 



dX £ t = v(X%)dt + y/edWt , I ( £ 6l 2 . (1) 

Here v(x) is an incompressible periodic vector field, W t is a 2-dimensional Brownian 
motion, and e (molecular diffusivity) is a small parameter. We further assume that the 
stream function if(xi,£ 2 ), such that 

V ± H=(-H' X2 ,H' Xi ) = v, 

is itself periodic in both variables, that is the integral of v over the periodicity cell is zero. 
For simplicity of notation assume that the period of H in each of the variables is equal to 
one. 

It is well known (see for example [4J), that with e fixed, the solution of ([1]) scales 
like a diffusion process with constant diffusion matrix when time goes to infinity. More 
precisely, there exists the limit, called the effective diffusivity, 

vei v £ j 

D«(e) = lim E A -^-^ , i,j = 1,2, 

t— >oo £ 



1 



where i and j are the coordinates and A is the initial distribution of the process Xf, 
which we can take to be an arbitrary measure with compact support. The measure 
on C([0, T], M 2 ), induced by the process -^X e ct ^ converges weakly, when c — > oo, to the 
measure induced by the diffusion process with constant matrix D(e). 

We are interested in the behavior of the effective diffusivity when the molecular dif- 
fusion e tends to zero. Assume that all the critical points of H are non degenerate. We 
distinguish two qualitatively different cases, depending on the structure of the stream 
lines of the flow given by v(x). 

In the first case, there is a level set of H, which contains some of the saddle points, and 
which forms a lattice in R 2 , thus dividing the plane into bounded sets, invariant under 
the flow. A standard example of a cellular flow, which has been studied in several of the 
papers cited below, is the flow with the stream function H(x\,X2) = sin(27rxi) sin(27TX2). 
In this particular example the separatrices (the level sets of H containing saddle points) 
form a rectangular lattice. 

In the second case, there is more than one unbounded level set of H containing critical 
points, and thus there are 'open channels' in the flow, and some of the solutions of the 
equation x'(t) = v(x(t)) go off to infinity. An example of a flow with open channels is the 
flow with the stream function H(xi, x 2 ) = sin(27rxi) sin(27rx 2 ) + 10 sin(27rx 2 ). Indeed, the 
horizontal axis {x 2 = 0} is an unbounded stream line of the flow. 




cellular flow flow with open channels 



Since v(x) is periodic, we may consider x'(t) = v(x(t)) as the flow on the torus. The 
torus is then a union of the sepatatrices and a finite number of open domains, bounded 
by the separatrices, and invariant under the flow. 

In [3] Fannjiang and Papanicolaou considered cellular flows for which the separatri- 
ces form a rectangular lattice on M 2 and the stream function satisfies certain symmetry 
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conditions. They showed that in this case 

D ij (e) = (d ij + o(l))y/e , as e^O, (2) 

that is the effective diffusivity is enhanced by a factor of order compared to case 
of the diffusion process \feW t without the advection term. Moreover, they found the 
constant matrix d %: > explicitly. Their proof is based on a variational principle applied to a 
symmetric operator associated to the generator of the process Xf. In j6] Heinze provided 
certain upper and lower estimates on the effective diffusivity in the case of cellular flows, 
for which the separatrices form a rectangular lattice on R . 

There are earlier physical papers (PQ, [8], [9]), arguing that the asymptotics in (j2J) is 
true for particular flows. Our first result is the rigorous proof of this statement for general 
cellular flows. 

Theorem 1.1. Assume that an infinitely smooth periodic stream function H(xi,X2) de- 
fines a cellular flow, and that its critical points are nondegenerate. Then the asymptotics 
of the effective diffusivity for the process (QP is given by 

Let C p be the noncompact connected level set of H. This level set contains some of the 
saddle points of H and forms a lattice in R 2 . Without loss of generality we may assume 
that H(x) = when x G C p . The corresponding level set on the torus will be denoted by 
C. 

The process X\ consists of the 'fast' part, which is the periodic motion along the 
streamlines, and the 'slow' diffusion across them. The motion is almost periodic away 
from the separatrices. However, once the trajectory is in a sufficiently small neighborhood 
of the level set £ p , it is likely to continue along it, and may go from cell to cell in a time 
much shorter than it would take the 'slow' diffusion to cover the same distance. 

The rough outline of the proof of Theorem II .11 is the following. We introduce a Markov 
chain, which can be viewed as a discrete time version of the process X\. The state space 
for the Markov chain is C. Note, that due to the periodicity of H, the process X\ can 
be viewed as a process on the torus. In order to define the transition probabilities, we 
introduce stopping times for the process Xf. The stopping time Tq is the first time when 
XI hits C, and is defined as the first time after r^_ x when the process X\ returns to 
C, after having traveled 'past' a saddle point. The transition times of the Markov chain 
are random. 

We show that the study of the asymptotics of the effective diffusivity can be reduced 
to the study of the asymptotics of transition probabilities and of the expectations of the 
transition times for the Markov chain. The limit of the transition probabilities as e — > 
is determined by the behavior of the process X\ in an arbitrarily small neighborhood 
of C The asymptotics of the expectations of the transition times, on the contrary, is 
determined by the event that the trajectory of X\ wanders away from the level set C 

In order to study the transition times we use the results of Freidlin and Wentzell 
[5]. For a given stream function H they introduce a graph and a mapping g from the 
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plane into the graph, such that each connected level curve of H gets mapped into a point 
on the graph, with the level sets containing critical points mapped into vertices. Then 
they demonstrate that the process g(X^ e ) on the graph converges to a limiting Markov 
process. The asymptotics of the expectations of the transition times for our Markov chain 
is related to the limiting process on the graph. 

Now consider the flows with 'open channels'. Assuming that the channels are directed 
along the x\ axis, we prove that 

D n (e) — (gP + j and D 22 (e) = (d 22 + o(l))e , as e — , (3) 

that is the diffusion across the channels is not qualitatively enhanced, compared with the 
process y/eWf The effective diffusivity in the direction of the flow is enhanced by a factor 
of order e~ 2 . 

Theorem 1.2. Assume that an infinitely smooth periodic stream function H{x\,X2) de- 
fines a flow with open channels, which are directed along the x\ axis, and that its critical 
points are nondegenerate. Then the asymptotics of the effective diffusivity for the process 
(QJj is given by (T3)]. 

Remark Since the matrix of effective diffusivity is symmetric and positive definite, from 
Theorem 1 1.2 1 it follows that the off-diagonal terms D 12 (e) = D 21 (e) are bounded uniformly 
in e. However we do not make a statement here on their asymptotic behavior. 

The proof of Theorem 11.21 is much simpler than that of Theorem 11.11 and is an easy 
application of the results of Freidlin and Wentzell [5] . 

The paper is organized as follows. In Section [2] we describe the construction of the 
discrete time Markov chain associated with a cellular flow, and relate the question of 
effective diffusivity for the process X f e to the study of transition probabilities and transition 
times for the Markov chain. In Sections [3] and H] we study of transition probabilities and 
transition times respectively. In Section [5] we prove Theorem 11.21 In Section [6] we prove 
several technical lemmas used in the previous sections. 

2 Construction of the Discrete Time Markov Chain 

Consider the set {x : < e ai } on the torus, with some j < a% < |, and let V s be 

the connected component of this set, which contains C Thus V s is a thin tube around 
C, whose width, however, is much larger than a typical fluctuation of the process H(X^) 
in fixed time (see the picture below). 

Let Ui, i — 1, n, be the connected components of T 2 \ £, and let A4, i — 1, n, be 
the saddle points, which belong to C. While the numbers of the connected components 
and of the saddles are the same for topological reasons, their equality is not used in the 
proofs, and the numbering of C/j's is not related in any manner to that of A^s. If there are 
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points, which are carried to Aj by the flow x' = v(x), and to Ai by the flow x' = —v(x), 
then the set of such points (a subset of C) is denoted by "f(Ai, Aj). We assume (for the 
sake of simplicity of notation only) that 7(Aj, Ai) is empty, that is the separatrices do not 
form "loops" . In a neighborhood of each curve 7(Aj, Aj) we may consider a smooth change 
of coordinates (xi,x 2 ) — > (H,9), where 9 is defined by the conditions: |V0| = \VH\ on 
7(Aj, ^), and W _L VH (this way 6 1 is defined up to multiplication by —1 and up to an 
additive constant). The same change of coordinates can be considered in V s D U k - In this 
case 9 G [0, f aUk \VH\dl], with the end points of the interval identified. Thus, if Ai G dU k , 
we define 

B(Ai, U k ) = {xeV £ nU k : 9(x) = 9{A i )}. 

LetB(A l ) = [J k . AtedUk B(A l ,U k ). 

Define the stopping times a £ = 0, Tq = inf{t : X s G £}. Then o"^,r^, n > 1 are 
defined inductively as follows. Assume that i G 7(A, Aj), and i 7^ j. Then 

< = inf{* > <„! : X t e G |J S(Afc) U^'}- 

Thus, o"^ is the first time after t^_ x that the process either exits V £ , or goes past a saddle 
point different from Ai. Define r £ = inf{t > o £ n : X £ G £}. Let C° = C \ {A^i = 1, ...,n}. 
Since almost every trajectory of X £ does not contain any of the points A iy X £ s is a Markov 
chain with the state space £°. The stopping times are the consecutive times when the 
process X\ hits the separatrix C after exiting V s or after having passed past a saddle 
point. Note, that the case, when a point x G ^{A-^Aj) travels, due to diffusion, against 
the flow v(x), and returns to B(Ai), does not count as having passed past a saddle point. 
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It is not difficult to see that X^ E satisfies the Doeblin condition, with the unique 
ergodic set. Therefore there exists a unique invariant measure (i £ (dy) on £° ([2]). Note 
that (X^ s ,t £ — Tn_i) also forms a Markov chain with the state space C° x R + , which 
satisfies Doeblin condition for each e, and has a unique ergodic set . Let p £ x (dy, dt) be the 
stochastic transition function for this chain (it only depends on the first component of 
the original point, as is reflected in the notation). Then Jl e (dy,dt) = J c0 p e x (dy,dt)fi e (dx) 
is the invariant measure. Since, due to the presence of diffusion, the distributions of the 
transition times have exponentially decreasing tails, we can apply the law of large numbers 
to the function (x, t) — > t on C° x R + to obtain 



lim = / / tjl £ (dy,dt) 

n-oo n J Jc xR+ 




tp^ x {dy 1 dt)^ £ {dx) = I E x T £ [i e (dx) almost surely. (4) 
IC°xC°xR + Jc° 

In the arguments, which led to (jl]), we considered X\ as a process on the torus. In order to 
keep track of the displacements of XL , as the process on R 2 , we introduce another Markov 
chain, on the extended phase space £ xZ 2 . Let now A{, % — 1, ... be the saddle points of H 
on the plane. Then any x p G 7(A> A?) c A> can ^ e uniquely identified with a pair (x, z), 
where x G C° and z = ([A}], [Af]) G Z 2 . ({A}] and [Af\ are the integer parts of the first and 
second coordinates of Ai). Thus, we have the mapping : C p \ {Ai,i = 1, ...} -4£°x Z 2 . 
Let 0i and 02 be the components of this mapping. We define the Markov chain Y £ as 
follows 

Note that the second component 02(X^ e ) — 2 (X^ e ) almost surely takes the values in 

some finite subset 5* of Z 2 . Thus Y£ is a Markov chain on C° x S, where S = {s\, Sk} C 
Z 2 . It is not difficult to see that Y£ satisfies Doeblin condition with the unique ergodic 
set. Applying the law of large numbers to the vector valued function f(x, s) = s, defined 
on C° x S, we obtain that there exists m G M 2 , such that 

X s e V" f(Y £ ) 

lim — — = lim — — — = m almost surely . 

n^oo n n^oo Ti 



Since — > N(0, D(e)) in distribution, and due to (HI), we conclude that m = 0. Applying 
the central limit theorem to the same function f(x, s), we obtain that there exists a matrix 
d e , such that 

X s V n f(Y £ ) 

lim — ^ = lim Z " i=1 i} 1 ' = N(0, d E ) in distribution . 



?woo y /n n—>oo x /n 

Due to (HI), we have 



D{e) = d e / / E x Tfd^(x) . (5) 
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In sections [3] and H] we shall obtain the asymptotics of the stochastic transition functions 
for the chain Y £ , and of the functions E x rf . Assuming that this is accomplished, the next 
lemma will allow us to obtain the asymptotics of the effective diffusivity, using ([5]). 

First we introduce the notations and formulate the assumptions needed for the lemma. 
Let M be a locally compact separable metric space. Let C b (M) be the set of bounded 
continuous functions on M. Let p e (x, dy), < e < e° be a family of stochastic transition 
functions on M. Assume that 

(A) The family of measures po(x, dy), x G K is tight for any compact set K. 

(B) p (x,dy) is weakly Feller, that is / f(y)p (x,dy) G C b (M) if / G C b (M). 

(C) For any / G C b (M) and any compact K C M, 

lim y f{y)p £ {x,dy) = J f(y)p (x, dy) uniformly in x G K . 

(D) There exist unique invariant measures /j, £ (dy). There exist A > 0, c > 0, such that 

\p n e (x, A) - fi £ (A)\ < ce" A " for all x,A,e . 

(That is p £ are uniformly exponentially mixing). 

Let g G C b (M, M 2 ) be such that J gdfx £ = for all e. Let Y £ be the stationary Markov 
chain, with the stochastic transition function p e . Since Y £ is exponentially mixing, the 
central limit theorem can be applied to g{Y £ ), and thus ^ ' converges weakly as 
n — > oo to a mean-zero Gaussian distribution. We denote the covariance matrix of the 
limiting distribution by d e (g). 

Lemma 2.1. Suppose that assumptions (A)-(D) hold. Then 
(a) // — > /i weakly. 

If f £ (x) are uniformly bounded, f°(x) G Cft(M), and lim^o / e (2 ; ) — uniformly 
on any compact, then 

lim y / £ d// = y . 

(c) If g G Cj,(M, M 2 ) zs snc/i t/iai J grf// = /or all e, then 

d £ (g)^d°(g) . 

Proof: From (A) it follows that for each n the family of measures p$(x,dy),x G K is 
tight for any compact set K. Let us assume that for a certain n, for any / G C b {M), 
uniformly on any compact set K we have 



f(y)p n e (x, dy) - J f{y)p n {x, dy) - as e . (6) 

Note that this is true for n = 1 by (C). Combining (jSJ) and the fact that Pq(x, dy),x G if 
is tight we obtain that for any compact set K and for any 5 > there is a compact set 
Ki such that 

p?(ar,tfi) > 1 — 5, x £ K (7) 
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for sufficiently small e. Next we justify ([6]) for n + 1 instead of n. 

f(y)(p £ l+1 (x,dy)-p n +1 (x,dy)) = 



f 1 (yM(x,dy)-p n (x,dy)) + J J 



f{z){Ps{y, dz) - po(y, dz))p"(x, dy) 



where fi(x) = f f(y)po(x,dy) G Cb(M). The first term on the right hand side tends to 
zero by (jSJ), while the second term tends to zero by (JTj) as f f(z)(p £ (y, dz) — po(y, dz)) is 
bounded and tends to zero uniformly on any compact. We therefore have established (jSj) 
for all n. 

Let us prove part (a) of the lemma. Fix an arbitrary x G M. Then for / G Cb{M) we 
have 

f(y)dfi e (y)= lim / f{y)p n £ {x,dy) 



and the limit is uniform in e by (D). The weak convergence of // to fj, now follows from 
©• 

To prove part (b) we write 

J fdfi £ - J f°dfi° = ( J fdff - J fdfi ) + ( J fdpf - J fdff) . 

The difference of the first two terms on the right hand side tends to zero as fi e —* fj, 
weakly. The difference of the last two terms tends to zero since for any 5 > there is a 
compact set K for which fJ, £ (K) > 1 — 5 for sufficiently small e (since // — > /i° weakly), 
and f £ — > f° on any compact set. 

In order to prove part (c) of the lemma it is sufficient to consider g G Cj,(M, R) (scalar 
valued). In this case 

d%g) = E[(g(Y £ )Y + 2g(Y*)g(Y{) + 2g(Y £ )g(Y £ ) + ...] . 

Due to uniform mixing (D) 

\E(g(Y £ )g(Y:))\ < , 

where 7 does not depend on e. In order to prove that d £ (g) — > d°(g) we therefore only 
need to establish that 

E[g(Y £ )g(Y:) - g(Y °)g(Y°)] - (8) 
for any fixed n. The left hand side of (JHD can be written as 

g(x)g(y)p n £ (x,dy)d^(x) - J g(x)g(y)p^(x,dy)d^(x) . 

Let G £ (x) = g(x) J g(y)p £ (x,dy). Then G e (x) — > G (x) uniformly on any compact by 
([5]), and the conclusion follows by part (b). □ 
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3 The Limit of the Transition Probabilities 



In this section we shall identify the limit of the transition probabilities for the chains Y£ 
on C° x S and verify the conditions (A)-(D) for these chains. 

Recall the (H, 9) coordinates which we may consider inside each cell Uk near its bound- 
ary, that is in Uk H V e . Let 

h = e-^H . 

In order to find the limit of the transition probabilities we shall demonstrate that in a small 
neighborhood of C p after a random change of time the process Xf is well approximated 
by the process X t with the generator \dhh + dg in (h, 9) coordinates. 

Let x G 7(Aj, A,) be a point on C p . We introduce the stochastic transition function 
Po(x, dy),x, y E C p as follows: 

Let x, y G dUk (if x and y do not belong to the boundary of the same cell, then 
Po(x, dy) = 0). Consider the (h, 9) coordinates in UkDV e , so that 6{x) = 0, and 9 increases 
in the direction of the flow. Since dUk is a closed contour, points with coordinates (h, 9) 
and (h,0 + J dUk \VH\dl) are identified. Let 9(Aj) and 9(y) belong to (0, f gu \VH\dl], and 
consider the process X t with the generator dg + \dhh in (h, 9) coordinates, which starts 
at the origin (corresponding to the point x). Let r be the time of the first exit from the 
following domain: V° = {9 < 9{A j )} U {9 > 9(A j ); h > 0}. Then define 

p (x,dy)= Prob a ,{fl(X T ) G [9{y)+m f \VH\dl,9(y+dy)+m f \VH\dl)}. 

(9) 

The summation over k is needed to account for the fact that x and y may both belong to 
the same edge j(A i: Aj), in which case they both belong to the boundaries of two cells, 
and we need to consider two sets of (h, 9) coordinates. 

The function po(x, dy) is a stochastic transition function on C p , and it can be consid- 
ered as a stochastic transition function on C° x S. It clearly satisfies conditions (A) and 
(B) preceding Lemma [2.11 

For x, y G C p , and for the stopping times defined in the previous section, let p £ (x, dy) 
be the transition function for the chain XL considered on C v : 

p £ (x, dy) = Prob x {X^ G [y,y + dy]} . 

Note that this definition is similar to that of po(x, dy). Here, however, we do not use the 
(h, 9) coordinates since with small probability the process X\ may travel outside of the 
domains Uk for which x G 7(v4j,A,) C dUk before time t[ (due to the presence of the 
small diffusion term X^ may go 'past' the saddle point Ai, thus traveling to one of the 
neighboring domains before time r-f). 

Lemma 3.1. For any closed interval I C "y(Ai, Aj), and any bounded continuous function 
f on C p , 

lim J f(y)p £ (x,dy) = J f(y)p (x,dy) uniformly in x G I . 
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Notice that Lemma 13.11 implies the condition (C) for the chain Y*. Before we start 
the proof of Lemma 13.11 we state and prove the following preliminary lemma. 

Lemma 3.2. Let X* and X 2 be the following two diffusion processes on M. d with infinitely 
smooth coefficients: 

dX\ = v(X})dt + a(Xl)dW t + e 2 v x {Xl)dt + ea^X^dWt , 

dX 2 = v(X 2 )dt + a(X 2 )dW t + e 2 v 2 (X 2 )dt + ea 2 (X 2 )dW t , 

with Xq = Xq. Suppose that for a certain constant L the following bound on the coeffi- 
cients holds: 

|Vi; < |,|Va iJ '|,|i;i|,|^|,|of|,|ai J '| <L , i,j = l,...,d, 

where i and j stand for the vector (matrix) entries of the coefficients. Let A be the initial 
distribution for the processes Xq and Xq. Then for some constant K = K{L) and for 
any t, rj > we have 

( e Kt _ i\ 2 

Prob A { sup |X] - X 2 \ > V } < { - . 

0<s<i V 

Proof: Let us assume that d = 1 in order to avoid vector and matrix indices. By Ito's 
formula, for any stopping time r < t, 

E x \Xl -X 2 \ 2 = E A [\(Xl - X 2 )(v(Xl) - v(X 2 ) + e 2 [ Vl (Xl) - v 2 (X 2 )])ds+ 
Jo 

E A f (a{Xl) - a(X 2 ) + e[ ai (Xl) - a 2 (X 2 )}) 2 ds . (10) 
Jo 

From the estimates on the coefficients and their derivatives it follows that the expression 
in (|TU|) can be estimated as follows 

E x \Xl - X 2 \ 2 < K(L)( [ E x \Xl -X 2 \ 2 ds + e 2 t). (11) 

Jo 

In particular, for r = t we have 



- X 2 \ 2 < K(L)( [ E x \Xl - X 2 \ 2 ds + e 2 t). 
Jo 



Let R(t) = E X \X} - X 2 \ 2 + e 2 . Then R(t) < K(L) J*R(s)ds, R(0) = e 2 . By Gronwall's 
Lemma applied to R(t) we have 

E x \Xl-X 2 \ 2 <e 2 (e K ^ t -l) . 

Define the stopping time r = min{s : \Xj — X 2 \ > rj} A t. Then, by (ITT|) 

r7 2 Prob A {max \\Xl - X 2 \ > r]} < E X \X^ - X 2 \ 2 < 

0<s<t 
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K(L)(E\ [ T \Xl - X 2 s \ 2 ds + eH) < e 2 {e K ^ - 1) , 
Jo 

which yields the lemma. □ 

Proof of Lemma I3.lt Since the kernel po(x,dy) is smooth in both variables and satis- 
fies condition (A) preceding Lemma [27T| to prove the uniform weak convergence stated in 
Lemma [3TT1 it is sufficient to demonstrate that for an arbitrary closed interval J C 7(A/, A{) 
and an arbitrary 5 > there is > such that 

p E (x, J) > po(x, J) — 5 for all x G /, e < Eq . (12) 

Suppose that I, J C dUk- (If / and J don't belong to the boundary of the same cell, 
then po(x, J) is equal to zero.) For the sake of simplicity of notation let us assume that 
/ C j(Ai,Aj) and J C j(Aj,Ai), that is I and J belong to the adjacent edges of C p . 
Without loss of generality we can assume that h > in £/& fl V e for sufficiently small 
e. We can consider the process X t e in (h, 9) coordinates in the following domain (see the 
picture below) 

V^ = {9> OiA,)} f|{>| < f| [{9 < 9{A ] )} \J{9 > 9(A,), h > 0}) . 




As above we consider the process X t in (h, 9) coordinates in the domain 

V° = {9< 0(Aj)} \J{9 > 6{Aj); h > 0} . 

Note that p e (x, J) is estimated from below by the probability that XI leaves T> £ through 
any of the copies of J (which corresponds to X\ making a finite number of rotations inside 
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Uk n V s , and then leaving Uk through the segment J). Let Jo, Ji, ... be the copies of J 
in (h,9) coordinates (J m +i can be obtained from J m by a shift by |Vif|cfZ along the 
9 axis). For an initial point x G T> e let p e (a?, J m ) be the probability that Xf leaves the 
domain T> £ through J m . As stated above, for x G / we have 



p £ (X 



J) > $^P £ (X, J m ) ■ (13) 



m=0 



Similarly, for a; G £>° let p"o(^> Jm) be the probability that X t leaves the domain T>° through 
J m . If x G I and J belong to different edges of dUk (as we have assumed) then by we 
have 

oo 

Po(aJ, •/) = y^Po(ff> Jm) • (14) 

m=0 

From the definition of po(x, J m ) it is easily seen that the sum converges uniformly in x 
for x & I. In order to prove (112p it is sufficient to demonstrate that for each m and each 
5 > there is e > such that 

p e (x, J m ) > po(x, J m ) - 5 for all x G J, £ < e • (15) 

Indeed, for every positive 5 we can select m such that po(x, J) < Sm=oP°( a '> Jm) + f • 
By ([15]) we can ensure that if we take Eq to be sufficiently small then 

Pa(x, Jm) < Pe(x, J m ) + 777 TTT for a11 X & 1,6 < £ , 171 < m , 

2(m + 1) 



which implies ( llzl) due to (113]) . 

For the sake of simplicity of notation we shall only prove (|15|) for m = 0. 
The generator of the process Xf is 

which in (h, 9) coordinates becomes 

L £ f = \WL\VH\ 2 + efnW\ 2 + y/if h AH + ef' 9 A9) + f e \VH\\V9\ . 

Dividing all of the coefficients of the generator by the same function |Vi/||W| amounts 
to a random time change for the process X%, which does not affect any of the transition 
probabilities. We shall denote the time-changed process with the generator L e f= ^jy^-, 
also by Xf. This process satisfies the equation 



dXl = (1, oW^|W + (0, l)V~eJ^dW e t + (4-^—, , + 1)* 



IV0I * v ' ' \ \VH\ 1 v 2 iVflllViZV 2 IV^IIV-H"! 



(16) 
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while 

dX t = (1, 0)dW t h + (0, l)dt , (17) 

where and Wf are one dimensional Wiener processes in h and 9 variables respectively. 
We can not apply Lemma [3.21 to (fTEj) and f[T7|) directly, as some of the coefficients on the 
right hand side of ( |T6l) may be unbounded near the saddle point (h = 0,9 = 9(Aj)). To 
circumvent this problem we shall take a sequence of steps (justified below), which will 
single out a small neighborhood of the saddle point, where yet another coordinate system 
will be considered. 

Step 1. Let us take 5' > small enough, so that there exist < hi < h 2 , such that 
for any 9\ G [9(Aj) — 5',9(Aj)] the probability of the event that X t passes through the 
interval K = {hi < h < h 2 ,9 = 9{\ before leaving the domain T>° through J differs from 
Po(x, Jo) by less than ^ for any x E I. Due to the smoothness of the transition kernel of 
the process X t , for some a < 1 the interval K can be replaced by any set K 1 as long as 
K' is contained in K and has Lebesgue measure at least (h 2 — hi)a. 

Step 2. If necessary make 5' from Step 1 smaller, so that |po(A Jo) ~Po(B, Jo) I < Jq 
whenever \h{A) - h(B)\ < 5' and A, B e [h u h 2 ] x [9{A 3 ) - 5\ 9(Aj) + 5']. 

Step 3. Take 5" < 5' and let K x = {h < h < h 2 , 9 = 9{A j )-5"}, and K 2 = {hi < h < 
h 2 ,9 = 9(Aj) + 5"}. Take 5" sufficiently small so that whenever A G Ki the process X\ 
starting at A passes through the set {h(A) - < h< h(A) + 9 = 9{A j ) + 5"} 

before leaving V s with probability at least 1 — ^ for small enough e. 

Step 4. Let us split the set K 2 into intervals Af , i = 1, ...,r of length 5' (we can 
assume that each of the intervals included the endpoints). Construct on Ki the intervals 
A} , i — 1, ...,r of length a8', such that ^(center (Aj)) = ^(center (A^)). Let tk x be the 
first time when a process either exits T> £ (T>°) or reaches K\. Let us take e§ so small that 

Prob^X? G A, 1 } > Prob x {X Tifi G A, 1 } - for all x G J, e < e . 
1 1 l(Jr 

Step 5. Let us take Eo so small that 

p £ (A, J ) > p (A, J ) - — for all A G K 2 , e < e . 

Assuming that the Steps 1-5 are valid let us prove (fl~5|) . By the Markov property 

r 

Pe{x,J ) > ^Prob s {X= G Aj} mm p £ {A, J ) . 

1 = 1 1 

By Steps 3,5, and 2, the second factor on the right hand side can be estimated as follows: 

S 2(5 3S 

mm p e (A, J ) > min p £ (A, J ) - — > min p (A, J ) - — > m&xp (A, J ) - — , 
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while by Step 4 

ProMX^ G A - } > Ymh x {X TKi G A - } - A 

Combining the above inequalities and using Step 1 we obtain 

8 w 3<> 

,(Frob4A rjfi G Aa-- 

8=1 



Pe(x, Jo) > ^(Prob^X^ G A, 1 } - — )(maxp (A Jo) ~ yjj) > 



VProb^X e Aj} maxp (A Jo) ~ ttJ > p(z, Jo) - 777 , 
i=i 

which implies (1151) . 

It remains to justify the construction in Steps 1-5. The validity of Steps 1 and 2 
follows from the fact that the transition kernel of the process X t is smooth. To justify 
Steps 4 and 5 it is sufficient to consider both processes X £ and X t in a compliment to a 
neighborhood of the saddle point, where Lemma [3.21 applies. 

In order to justify Step 3 we note that by Morse Lemma in a neighborhood Oj of 
the saddle point Aj there is a smooth change of variables, such that in the new variables 
the stream function is H{xi,x%) = X1X2, and the interior of Uk corresponds to the first 
quadrant Xi, X2 > 0. In the new variables the generator of the process X^, after a random 
change of time, becomes L £ f = eLif + v{Vf, where L\ is a differential operator with first 
and second order terms, with bounded coefficients, and V\{x\,x-i) = (— X\,Xz). We shall 
consider the operator L £ in the domain V s = Oj f]{x\ > 0; X2 > 0; x\ + X2 > ; X1X2 < 
e^}. Make a further change of variables in T> £ : 

(xi, x 2 ) -> (u, v) = {^-p-, x 2 - Xi) . 

In the new variables, after dividing all the coefficients of the operator by {x\ + X2), which 
amounts to a random change of time for the process, the operator can be written as 

Uf = M £ f + ^- ) (18) 

where M e f is a differential operator with first and second order terms. All the coefficients 
of M £ can be made arbitrarily small in T> £ by selecting a sufficiently small neighborhood 
Oj of the point Aj, and then taking e to be sufficiently small. 

The construction in Step 3 now follows from Lemma 13.21 by comparing the process 
whose generator is the operator ( fl8i) with the deterministic process with generator |£. 

Remark To verify condition (D) (uniform mixing) preceding Lemma 12.11 for the chain 
Y £ it is sufficient to show (see [2], page 197) that there is an integer n > 1, an interval 
I G C° x S, and a constant c > 0, such that 

p"(x,dy) > c\(dy), and Pq(x, dy) > c\(dy), for x G £° x S, yel , (19) 
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where X(dy) is the Lebesgue measure on I. The proof of estimate (fT9l) is absolutely similar 
to that of Lemma 13.11 



4 The Asymptotics of the Transition Times 



In this section we shall study the asymptotics of the integral J £0 E^rf <i/i e (x), which enters 
in the expression ([5]) for the effective diffusivity. 



We shall demonstrate the following: 




lim^E^Tf = f°(x) uniformly in x G /, (21) 



where f°(x) G C\,{C°) is a positive function and I is an arbitrary closed interval / C 
j(Ai,Aj). From parts (a) and (b) of Lemma [2. II it then follows that 



where fi°(x) is the invariant measure on £° for the kernel po(x, dy), defined in Section [31 
The proof of formulas fl20l and ff2T]) will rely on a sequence of lemmas stated below. 
We shall study separately the probability of the event that the process X| starting form 
x G j(Ai, Aj) reaches dV e before time r-f , and the expectation of the time it takes for the 
process starting from dV e to reach CP. 

Consider the process X\ together with the process X t , whose generator in (h, 9) co- 
ordinates is \d hh + d e in the domain T>\ = {9{Ai) < 6 < 9(Aj); \h\ < e ai '^}. We follow 
the process X t till it exits T>\ . Let ^0(^5 dh) be the corresponding transition kernel. Thus 
Po(x,dh) coincides with a Gaussian distribution on — e ai ~2 < h < £ ai_ 2 ; and has two 
point masses at h — ±e ai ~2. Similarly let P £ (x,dh) be the transition kernel for the 
process Xf, which starts at x G j(Ai, Aj) and is stopped at the time a\. We have the 
following: 

Lemma 4.1. For any continuous function / : R — > R, such that \f{h)\ < 1 + \h\, there 
exists c> 0, such that 




E x Tf dfi e (x) =£~H[ f°{x)dfj,°{x) + o(l)) as e -> , 



Jc° 



(22) 




(23) 



Furthermore, for any closed interval I C j(Ai, Aj), 




uniformly in x G I . 



(24) 
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The proof of Lemma 14.11 is completely similar to that of Lemma 13.11 
We introduce the following notation: T VE is the first time the process XI leaves V £ ; 
similarly, r Uk and the first instances when X\ leaves Uk and V s fl Uk respec- 

tively. 

In order to estimate the probability that the process Xf starting from x G 7(^4,, Aj) 
reaches dV £ before time rf we shall need the following 

Lemma 4.2. FJ]/ There exists a constant c > 0, such that 

E x r y£ < ce 2ai ~ l \\ne\ for any x G V s . 

This Lemma is the same as Lemma 4.7 of [5] (it must be observed that the proof of 
Lemma 4.7 of [5] goes through for any A24 < |). 

In the event that X £ e G V s fl Uk, after the stopping time <rf the process Xf may exit 
V s fl Uk either through dV £ or through C°. The next lemma estimates the probability 
that the process exits the domain through dV £ . 

Lemma 4.3. There exists c > 0, such that for any x G V £ C\Uk 

\Prob x {X £ v e nUk G dV £ } - h(x)e^~ ai \ < ce ai \ lne\ . (25) 

Proof: Let L £ be the generator of the process X £ in the domain V £ C\ Uk- Then the 
probability in (TS5]) is equal to the solution u(x) of the equation Uu = in V £ fl Uk with 
the boundary conditions u\H= £ a i = 1, u\h=o = 0. Let U\{x) = u{x) — -f^. Then u± is the 
solution of the equation Ifu\ = —L 6 ^^- with the boundary conditions ^i|a(y £ nc/ fc ) — 0. 
By Lemma [4.2[ since e~ l L £ H(x) is bounded uniformly in e, the solution u\ is estimated 
as follows: 

|«i| < c Q e l - ai V x T™ k < Cl e ai \\ne\ . 

This implies the statement of the lemma. □ 
Using the Markov property of the process X £ with respect to the stopping time a{, 
we obtain that there is c > 0, such that for any x G C° we have the following: 

/oo 
sup _Prob,{X r %. en ^ ) G dV £ }P £ (x,dh) < ce^ , (26) 
■00 x&V £ :h(x)=h 

where U(x) is the domain which contains x (one of the domains Uk), and the second 
inequality is due to Lemma fl~3l and ( 1231) . Furthermore, due to ( 1241) we can evaluate the 
asymptotics of the event {r ye < rf ; X E ve G dV £ fl Uk} as follows: 

POO 

limProb x {r y£ < t £ ; X £ vs G dV £ T)U k }e ai ~^ = lim / hP (x,dh) , uniformly in x £ I. 

(27) 

The next lemma allows us to estimate the expectation of the time it takes for the process 
starting at dV £ to return to C°. 
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Lemma 4.4. For each of the domains Uk there exists a constant Ck > 0, such that 

\ime 1 ~ ai 'E x r Uk =c k uniformly in xedV £ f)U k . (28) 

Note that Lemma 14.41 together with (1261) and (1271) implies (1201) and (12ip since the 
expectation of the time it takes for the process to reach dV £ can be estimated by Lemma 
14.21 It remains to prove Lemma [4.41 

We introduce notations and state several technical lemmas needed for the proof of 
Lemma 14.41 

Since H(x) = on dUk, we may assume without loss of generality that H(x) > 
inside Uk in a small neighborhood of Uk- Then there is a region V C Uk, whose boundary 
consists of dUk and a level curve {H(x) = H }, and, by selecting a sufficiently small H , 
we can ensure that each level set of H in V is connected and there are no critical points 
of H in the closure of V other than on dUk- 

For < H < Ho let us define the following functions: 

a(H) = J \VH\dl, b(H) = J -^-dl, q(H) = J — Ldi , (29) 

in each case the integration is over the level set {H(x) = H,x £ V}. Let r < ^ be a small 
constant, to be specified later. Consider the function f(H), which solves the equation 

a(H)f(H) + b(H)f(H) = -q{H) , (30) 

with boundary conditions /(0) = /(2r) = 0. While it not used here explicitly, we note 
the fact that the operator in the left hand side of (!30|) after dividing it by the function 
2 q{H) becomes the generator of the limiting diffusion process on the edge of the graph 
corresponding to the domain Uk (cf [S] and Section Oof this article). 
We need the following lemma, which will be proved in Section 

Lemma 4.5. There is a function g(r) , which satisfies lim r ^ g( r ) = 0, such that \f'(H)\ < 
g(r) for all < H < 2r. Further, there is a constant c > such that \f"(H)\ < c\ lnH\ 
and\f"(H)\<§. 

Let us select constants a 2 and a 3 such that a\ < < a 3 < |. Define the subsets V A 
and V B of V as follows: 

V A = { x e V; e a2 < H(x) < r}, V B = {x G V; e az < H{x) < 2r} . 

Let t a be the first time the process X s leaves V A , similarly r B is the first time the process 
leaves V B . Let xt be the deterministic process 

dx t = v(xt)dt , 

and let T(x) be the time it takes the process x t starting at x to make one rotation along 
the level set, T(x) = mf t>0 {x t = x}. The next lemma shows that for times of order T(x) 
the process XI is in a certain sense close to the deterministic process Xt- The lemma is 
proved in Section El 
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Lemma 4.6. For any 5 > there is 7 > such that 

Prob x { sup \H(X £ S ) - H(x s )\ > e^ 5 } < for all x G V" 4 . (31) 

s<T(x) 

There exist 5' > and 7 > suc/i £/ia£ 

Prob^ sup \X £ S - x 8 \ > e s '} < e 1 for all x G V A . (32) 

s<T(x) 

One of the main ingredients of the proof of Lemma 14.41 is the following lemma, which 
is a particular case of the main result (Theorem 2.3) of [5]. 

Lemma 4.7. |5j/ There is a constant > 0, such that on each level set {H(x) = r, x G V} 
we have 

limeE^* =c fc (l + (/i(r)) , (33) 

the limit is uniform on each level set, and gi(r) satisfies lim r ^ 9i{ r ) — 0. 

Lemma [4.71 is different from Lemma [4.41 in that the initial point in (1331) belongs to a 
fixed level set of H, while in (128]) the initial point is asymptotically close to dU^ as e — > 0. 

Proof of Lemma 14. 4t We shall demonstrate that there exists a function g(r), such 
that lim r ^o5'( r ) — 0, for which 

E,r A < e a ^ l g{r) , (34) 

uniformly in x G {H(x) = e ai ,x G V}. Let us show that (|34"1) is sufficient to prove the 
lemma. As in the proof of Lemma 14.31 from ( 1341) it follows that 

iProb^if^) = r} - — | < c— <?(r) (35) 

uniformly in 2 G {-ff(ar) = £ ai , x G V}. For x G <9V £ H by the Markov property 

E x r u * = E x r A + E x (E xu t u "; H(X* A ) = r) + E^ (E^ A r Uk ;H {X^ A ) = e a *) . (36) 

The first term on the right side of (1361) is estimated from above by e ai ~ 1 g(r) due to (1341) . 
The second term has the following asymptotics due to (131)1) and Lemma 14.71 

lE^E^r^; H(X* A ) = r) - c^^ 1 ! < ^(r)^ 1 " 1 

T 

for sufficiently small e, where g^O") satisfies lim r _»o fl^C?") = 0. The last term on the right 
side of ( 136]) is estimated from above by ce a2 ~ l due to Lemmas 14. 21 and 1 4.31 an d the repeated 
use of ( |36|) . Therefore 

\E x t u » - c k e a >~ l \ <e a -- l gz{r) 
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for sufficiently small e, and g${r) which satisfies lim r ^o 9z{ r ) — 0. Since r can be selected 
arbitrarily small we obtain lim e ^ = c&. It remains to prove (13~4"1) . 

Let a B = min{r B , T(x)}. We shall prove that for some K > for all sufficiently small 
values of r 

+ KVJWXU)) < Kfmx) ) 

e e 

for all x G V A . From ( 1371) it follows that 

E x r A < Kf{ f x)) for x G V A . (38) 

Due to the estimate on the derivative of / from Lemma [4.51 for x G {H(x) = e ai ,x G V} 
the right side of (138]) is estimated from above by Kg^e 011 ' 1 , which implies (l34"j) . Now we 
need to prove (l37j) . 

Applying Ito's formula to f(H(Xf)) we obtain 

V/(tf(X^))-/(tf(*))) = ~E S ^ (f"(H(X!))\VH(X!)\"+f'(H(Xl))AH(X!))ds, 
£ ^ Jo 

while from (130]) 

T(x) = - [ T(X \f"(H(x s ))\VH(x s )\ 2 + f(H(x s ))AH(x s ))ds . 
Jo 

Thus, what we want to show is that there is a constant K\ such that for all x G V 



A 



E x \ ^ [f"(H(Xt))\VH(X^)\ 2 - f"(H(x s ))\VH(x s )\ 2 ]ds\ + 
Jo 

E x | / [f'(H(X!))AH(X!)) - f'(H(x s ))AH(x s ))]ds\+ (39) 
Jo 

E,| / TW [r(^(Xf))|V J ff(Xf)| 2 + f(^(^))A/f(X|))]^| <K x T(x) . 

Since /' and Aif are bounded, and K\ can be taken arbitrarily large, it is sufficient to 
estimate only those of the terms in (|39"1) which contain the second derivative of /. By 
Lemma 14.61 we have Prob^cr 5 < T(x)} < e 7 , and the second derivative of / can be 
estimated by Lemma 14.51 Therefore for the last term containing /" we have 

x \ [ T(X) f\H(Xl))\WH(Xl)\ 2 ds\ < cT(x)\ \n(e a3 )\Pvob x {a B < T(x)} < cT(x) . 

Ja B 



E 

The estimate 



E x | f [f"(H(X!))\VH(X!)\ 2 - f{H{x s ))\VH{x s )\ 2 ]ds\<cT{x) 
Jo 



19 



follows from Lemma 14.61 and the estimates on /" and /"' of Lemma 14.51 This completes 
the proof of Lemma 14.41 □ 



Proof of Theorem ll.lt The effective diffusivity D(e) is related to the variance d £ 
of the limit of the functional of the Markov chain, and to the integral of the expectation 
of the transition times via formula (j3J) . As shown in Lemma 13.11 and in the Remark fol- 
lowing it, Lemma 12.11 applies, and therefore there exists the limit d° = lim e ^o^ £ - The 
asymptotics of the integral f~ E x r[dfi e (x) is given by fl22|) . This completes the proof of 
the theorem. □ 



5 The Case of the Open Channels 

In this section it will be convenient to consider the process X^, which is the same as X%, 
but only accelerated by the factor |, that is Xf = X^ £ . This process satisfies the equation 

dX £ t =-v (XI) dt + dW t , X s t e R 2 . 

e 

Note, that as a process on the torus, Xf is uniformly (in e) exponentially mixing. Fol- 
lowing [5J we consider the finite graph G which corresponds to the structure of the level 
sets of H on the torus. 

The graph G is constructed as follows: we identify all the points which belong to each 
connected component of each level set of H. This way each of the domains Uk, bounded by 
the separatrices, gets mapped into an edge of the graph, while the separatrices themselves 
get mapped into the vertices. Let e(Xf) label the edge of the graph and let if(X t £ ) be 
the coordinate on the edge. Then the process (e(X^), H(Xf)) can be considered as a 
process on the graph. It is proved in [5] (Theorem 2.2) that the process (e(Xf), H(Xf)) 
converges to a certain Markov process on the graph with continuous trajectories, which 
is exponentially mixing. We state the result here in less generality than in [5], but this is 
sufficient for our purposes. 

Theorem 5.1. /3J/ There is a Markov (diffusion) process Y t on the graph G, which is 
exponentially mixing, and has continuous trajectories, such that for anyT > 0, the process 
(e(Xf), H(Xf)) converges to Y t weakly in C([0,T],G). 

We are now in the position to prove Theorem 11.21 
Proof of Theorem 11.21 The displacement of the process X\ in the direction x\ (the 
direction of the channels) is given by 

{Xlf = - f Vl {Xl)ds + Wl 
^ Jo 
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where v\ is the %\ component of the velocity field. Therefore, 



(i) = ite « +o(1)) 



Exv^Xfiv^X^ds + o(l)) , 



where Xq is distributed according to the invariant (Lebesgue) measure A on T 2 . For a 
function / G C°°(T 2 ), let /(e, if), (e, if) G G be the function defined on the graph, other 
than on the vertices, which is equal to the average of / over the corresponding connected 
component of the level set of H 

lie H) - ^ fiXs)dS 

where x t is the solution of the deterministic equation dx t = v(x t )dt, the initial point x 
belongs to the level set, and T(x) is the time of one revolution around the level set. It is 
easily seen that for any initial point x which does not belong to any of the separatrices 
of H we have 

lim / E x f(X!)ds = 0, if J(e,H) = 0. 

Therefore, 

lim [ t E x v 1 (XZ)v 1 (XZ)d8 = Vm f E^X^MK), H(Xf)) ds = 
Jo Jo 

" E.v^v^ds , (40) 



I 

Jo 



'0 

where \i is the measure on G, which is invariant for the process Y t . 
The integrals 

/oo 
E xVl (X £ )vi(X £ s ) ds 

and 

/oo 
E^Yop^Ys) ds 

can be made arbitrarily small by selecting sufficiently large t due to uniform mixing of 
the processes Xf and Y t . Therefore, 



lim / E x v 1 (XZ)v 1 (X e a )ds= / E li v 1 (Y )v 1 (Y 8 )ds, (41) 
e ^ u Jo Jo 

which shows that the asymptotics for D u (e) is as stated in the theorem. 
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Now let us consider the asympotics for D 22 (e). Note that ^(e, H) = 0, thus the 
arguments leading to (HTj) do not provide the asymptotics of D 22 (e). Let P 1; . . . ,P n be 
those of the separatrices of H on the torus which, when unfolded onto the plane, are 
non-compact. Let us select a point Ai on each of Pi. Let us introduce the sequence of 
stopping times r n , n > 1, which are the consecutive times when X s makes the transition 
to a different level set Pi. Thus X e T is a Markov chain on the set {Pi, . . . , P n }. We can 
also consider the Markov chain 

%n = (X!; n ,T n — r n _!, A n ) 

on the extended phase space {Pi, . . . P„} x M + x R. The third component A n is defined 
as follows: If XI is considered on the plane, then 

A n = A 2 (n)-A 2 (n-l), 

where A 2 (n) is the x<i coordinate of the point corresponding to the separatrix containing 
the point X £ n . 

Similarly we can introduce the stopping times rj n for the process Y t on the graph, which 
are the consecutive times when Y t visits different vertices Qi = H(Pi) of G, corresponding 
to the unbounded separatrices of H. Together with the Markov chain Y Vn we can consider 
the chain 

on {Qi, . . . Q n } x M + x R, where A n is defined the same way as A n . 

Let fi £ be the invariant measure for the chain Z e nl and let jl be the invariant measure 
for the chain Z n . Let / be the function defined on the state space of the chain Z e w which 
is equal to the third component: f(x, r, A) = A. The function / is defined the same way 
on {Q u ...Q n } x K + x R. 

By the central limit theorem applied to the chain Z £ w there is a number d e such that 



Similarly, 



lim YH=of{Zi) =N , Q d ) m 

n^oo n 



The effective diffusivity in the %2 direction is then different form d e by the factor 



,22,^ £d ( £ ) 



From Theorem 2.2 in [5] and the uniform mixing of the Markov chains Z £ n and Z n it 
easily follows that 

d(e) — > d 
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and 

T\d[i £ — > / 77i<i/I. 



This completes the proof of Theorem 11.21 □ 

6 Proof of the Technical Lemmas 

Proof of Lemma 14.51 For any function u G C°°(Uk) we have 

dHj udl = J ww J — ivffi — (42) 

where the integrals are over the level set {H(x) = H,x G V}. In particular b(H) = a'(H), 
and therefore equation (130]) can be written as 



(a(H)f(H))' = -q(H) . (43) 

From the definition of the coefficients a(H), b(H), and q{H) it easily follows that 

lim a(H) = a > ; b{H) = 0(| ln#|) as H ^ ; q(H) = 0(\ \nH\) as // -f . 

(44) 



Further, with the help of Morse Lemma and (142]) it is easily seen that 

b'{H) = O(^r) as H -> ; g'(iJ) = 0(4) as # -> . (45) 

Let H m G (0, 2r) be the point where f(H) achieves its maximum, thus f'(H m ) = 0. From 
( 143]) it follows that 

— fu q(s)ds 

f(H) = JH f\ J . (46) 

Thus, the estimate on the first derivative of / stated in the Lemma follows from (jHJ). 
Rewrite (130]) as 

q(H) + b(H)f'(H) 

f {H) = ajH) ' 

From (jUJ) in now follows that \f"(H)\ < c|ln_ff| for some c > 0. Differentiating both 



sides of ( 130]) we obtain 

q'(H) + b'(H)f(H) + b(H)f'(H) + a'{H)f"{H) 



f"{H) 



a(H) 



The estimate on f"'(H) now follows from the estimates on the first two derivatives and 
from (T441 and (143]) . This completes the proof of the Lemma. □ 
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Proof of Lemma 14.61 The proof is based on the use of Lemma 13. 21 We can not 
however apply Lemma 13.21 to the pair of processes Xf and x t directly, since the rotation 
time T{x) grows logarithmically in e when x G V A . 

Let us establish the following property of the deterministic flow xf. 

Let = t < ti < t 2 < ... < t n . Consider a process y t , which solves the equation 

dy t = v(y t )dt (47) 

on each of the segments [t ,ti), [ti,t 2 ), [t n -i,t n ], with a finite number of jump discon- 
tinuities lim 4 ^ t . + y(t) — lim^f^ y(t) = Pi, i = 1, ...,n — 1. Then for any positive c there 
are positive k and 5' such that under the conditions 

x t0 = yto G V A ; ^2\\pi\\ < e^~ K ; t n - 1 < c\\ne\ 

i=l 

we have 

sup \\y t -x t \\< e 25 ' . (48) 
o<t<t„ 

Note that it is sufficient to establish the following: for any pair of points ao, b such that 
a G V A and ||a — &o|| < we have 

sup \\a t - b t \\ < E 28 ' +K ~^\\a - 6 || , (49) 

0<*<c| lne| 

where a t and b t are the solutions for the deterministic flow (JITI) . Let us take 

8' = k=\{\-<* 2 ). (50) 



The time it takes for the trajectory of fl47|) to make one rotation along the level set 
{H(x) = H, x G V} is equal to T(x) = q(H(x)), where q(H) was defined in fl29|) and is a 
smooth function for sufficiently small positive H, which satisfies 

q(H) = 0(| \nH\) , q\H) = 0(1) as H -+ . (51) 

The number of full rotations of the trajectory starting from ao is equal to [ y AJ and the 

time it takes to make these rotations is equal to [ r /J T(ao). It takes [ T ^ ]T(bo) to make 
the same number of rotations for the trajectory starting at bo. 
Due to ( IBTl) the difference is estimated as follows; 

|[^]T(a )-[-l T ]T(6o)l <oonst|lne|e-°»||ao-6d|| • (52) 
T(a ) T(a ) 

Here we used the facts that t < c\ In e| and that H(ao) > £ Q2 - Now consider the images 
of ao and bo under the flow ( 14T1) for time t < T(a ). Using the reduction of the flow to a 
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linear system in a neighborhood of each of the saddle points (Hartman-Grobman Theorem 
[7]), it is easy to show that 

sup \\a t — b t \ \ < const||a — &olk~ a2 • 

0<T(a ) 

Combining this with (|52j) and with the fact that the speed of motion in (|47p is bounded, 
we obtain (|49p with 5' and k defined in ( 1501) . This in turn implies (|48|) as noted above. 

Note that < 25' < | — k and that by making k smaller (if necessary) we can satisfy 
< k < 5, where 8 is the same as in (EH). Observe that for some c > we have 



T(x)<c\\ne\ for all x G V A . (53) 

Select the points = to < t\ < t 2 < ... < t n = T(x) in such a way that ^|ln£| < 
\ti + i — ti\ < jf \ In e| for i = 0, n — 1, where K is the constant from Lemma T3.2I (applied 
to the pair of processes Xf and Xt)- By (153]) there is the estimate n < Let yf 

be the piecewise continuous process, which is defined by the conditions: yf_ = Xf. and 
dy\ = v(yf)dt on [tj, t i+ i), i = 0, n — 1. By Lemma I3~21 

Prob x {V sup \\X £ t - y E t \\ > e^ K } < ( )V . (54) 

Due to continuity of Xf formula ( |54l) provides an estimate on the sum of the jumps of 
the process y\ . From (08]) it now follows that 

IrK 

Prob^Hxt-yfll >e 2S '} < ( )V . (55) 

K 

This, together with (154"]) implies (132]) for any 7 < k. Since H(x t ) is constant and H{yf) is 
piecewise constant, we have 

Prob^ sup \H(Xf)-H(x t )\>e*- s }<Prob x {J2 sup ||X f e - y f e || > } . 

t<T(x) ^mu,U+i) SUp ||Vi^|| 



This, together with (]54|) and the condition that k < 5 implies ( 1311) for any 7 < k. This 
completes the proof of the Lemma. □ 

Acknowledgements: I am grateful to Prof. S. Molchanov for introducing me to 
this problem and for many useful discussions. 
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